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Abstract  The  dynamics  of  fluid  particles  on  cylindrical  manifolds  is  investi¬ 
gated.  The  velocity  field  is  obtained  by  generalizing  the  isotropic  Kraichnan 
ensemble,  and  is  therefore  Gaussian  and  decorrelated  in  time.  The  degree  of 
compressibility  is  such  that  when  the  radius  of  the  cylinder  tends  to  infinity 
the  fluid  particles  separate  in  an  explosive  way.  Nevertheless,  when  the  radius 
is  finite  the  transition  probability  of  the  two-particle  separation  converges  to 
an  invariant  measure.  This  behavior  is  due  to  the  large-scale  compressibility 
generated  by  the  compactification  of  one  dimension  of  the  space. 
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1  Introduction 

Many  physical  systems  display  a  strong  dependence  on  the  space  dimension¬ 
ality,  the  best  known  example  being  given  by  phase  transitions  in  equilibrium 
statistical  physics.  As  for  non-equilibrium  systems,  hydrodynamic  turbulence 
shows  a  remarkable  dependence  on  the  space  dimension  as  well.  In  three  di¬ 
mensions,  the  kinetic  energy  flows  from  large  to  small  scales  in  the  form  of 
a  Kolmogorov  Richardson  cascade.  Conversely,  in  two  dimensions,  energy  is 
transferred  upscale  at  a  constant  rate,  in  an  inverse  cascade  process  |18j.  Ad¬ 
ditionally,  three-dimensional  turbulence  is  characterized  by  a  breakdown  of 
scaling  invariance  and  small-scale  intermittency  m ,  whereas  the  inverse  cas¬ 
cade  is  apparently  self-similar  [2]  and  even  shows  some  intriguing  signatures 
of  conformal  invariance  [1]. 

These  observations  have  spurred  the  search  of  a  critical  dimension  be¬ 
tween  d  =  2  and  d  =  3  in  the  hope  that  it  could  provide  a  starting  point  for 
an  analytical  attack  of  three-  or  two-dimensional  turbulence,  or  both.  This 
approach  has  been  mainly  applied  to  simplified  models  of  turbulence,  such 
as  EDQNM  [1 2].  the  shell  model  IB),  or  a  model  obtained  by  generalizing 
the  form  of  the  two-dimensional  stream  function  [24] .  In  those  studies,  the 
spatial  dimension  has  been  most  conveniently  reduced  to  a  formal  parameter 
that  could  take  arbitrary  values.  The  approach  undertaken  in  the  present 
paper  differs  from  previous  work  at  least  in  two  important  aspects.  First, 
we  shall  consider  a  geometrical,  rather  than  formal,  way  of  looking  in  be¬ 
tween  dimensions.  Namely,  we  shall  study  the  dynamics  of  fluid  particles  on 
cylindrical  manifolds  where  the  compactified  dimension  can  be  collapsed  or 
inflated  at  will  so  as  to  connect  continuously  the  two  extreme  cases.  Second, 
we  shall  focus  on  a  system  that  is  fully  under  analytical  control,  that  is  the 
Kraichnan  ensemble  of  velocities  rather  than  Navier-Stokes  turbulence. 

To  study  the  turbulent  transport  of  a  passive  scalar,  Kraichnan  intro¬ 
duced  a  Gaussian  ensemble  of  decorrelated-in-time  velocity  fields  M-  A 
compressible  generalization  of  the  Kraichnan  ensemble  in  the  d-dimensional 
Euclidean  space  has  been  investigated  by  Gaw§dzki  and  Vergassola  under 
the  assumption  of  statistical  isotropy  |T5]  0  In  this  model,  the  dynamics  of 
fluid  particles  depends  on  three  physical  quantities:  the  space  dimension,  the 
degree  of  compressibility,  and  the  (spatial)  Holder  exponent  of  the  veloc¬ 
ity.  The  Holder  exponent  £/2  is  greater  than  zero  and  less  than  one.  This 
property  mimics  the  behavior  of  a  turbulent  velocity  field,  whose  realizations 
are  typically  non-Lipschitz  in  the  limit  of  infinite  Reynolds  number.  For  any 
given  d  <  4  and  0  <  £  <2,  Gaw§dzki  and  Vergassola  have  identified  a  critical 
degree  of  compressibility  separating  two  different  phases  of  the  Lagrangian 
dynamics.  Below  the  critical  value  (incompressible  or  weakly  compressible 
velocity  fields),  fluid  particles  separate  superdiffusively.  The  probability  dis¬ 
tribution  of  fluid-particle  separations  does  not  have  a  stationary  limit  in 
this  regime.  Above  the  threshold  (strongly  compressible  fields),  Lagrangian 
trajectories  tend  to  collapse  to  zero  distance,  and  the  distribution  of  the  sep¬ 
arations  degenerates  into  a  Dirac  delta  function.  For  d  >  4,  the  former  regime 


1  The  smooth  limit  of  this  model  had  been  previously  considered  in  ref.  [3]. 
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is  the  only  possible  one  and  the  phase  transition  does  not  occur 0  The  above 
results  have  been  subsequently  elaborated  by  Le  Jan  and  Raimond  in  the 
context  of  non-Lipschitz  stochastic  differential  equations  mm- 

Here  we  consider  a  generalization  of  the  Kraichnan  ensemble  on  a  cylin¬ 
drical  surface.  A  d-dimensional  cylindrical  surface  can  be  constructed  by 
taking  Rd  and  compactifying  d  —  d!  dimensions.  The  radius  of  the  cylinder 
is  the  size  of  the  compactified  dimensions.  When  the  radius  tends  to  infinity 
we  recover  Rd;  when  it  tends  to  zero  we  obtain  Rd  .  Thus,  varying  the  ra¬ 
dius  of  the  cylinder  produces  a  smooth  transformation  from  dimension  d  to 
dimension  d! . 

We  define  a  zero-mean  Gaussian  velocity  field  on  a  cylindrical  surface  by 
imposing  the  form  of  its  covariance.  We  require  that  the  covariance  of  the 
field  tends  to  the  one  of  the  isotropic  d-dimensional  Kraichnan  ensemble  as 
the  radius  of  the  cylinder  tends  to  infinity  and  to  the  one  of  the  isotropic 
d'-dimensional  Kraichnan  ensemble  as  the  radius  vanishes.  The  degree  of 
compressibility  is  such  that  the  velocity  is  weakly  compressible  in  the  limit  of 
infinite  radius  and  strongly  compressible  in  the  opposite  limit.  It  is  therefore 
possible  to  gradually  move  from  one  regime  to  the  other  by  varying  the  radius 
of  the  cylinder. 

As  we  shall  see,  if  in  the  limit  of  infinite  radius  the  Holder  exponent  is 
equal  to  £,  then  in  the  limit  of  vanishing  radius  it  is  equal  to  £  +  (d  —  d'). 
Hence,  if  attention  is  restricted  to  non-smooth  velocities,  the  model  under 
consideration  is  meaningful  only  when  a  single  dimension  is  compactified 
(d'  =  d  —  1).  For  the  sake  of  simplicity,  we  shall  conduct  the  analysis  in  two 
dimensions,  where  the  two  extreme  cases  are  the  two-dimensional  plane  and 
the  straight  line.  We  shall  show  that  the  dynamics  of  fluid  particles  results 
from  two  opposite  effects.  At  small  separations,  Lagrangian  trajectories  ex¬ 
hibit  a  superdiffusive  dynamics  owing  to  the  weakly  compressible  nature  of 
the  small-scale  velocity.  At  large  separations,  fluid  particles  experience  the 
trapping  effect  of  a  strongly  compressible  field.  Consequently,  the  probability 
distribution  of  the  two-particle  separation  tends  to  an  invariant  measure.  This 
behavior  is  to  be  contrasted  with  the  one  observed  in  the  two-dimensional 
isotropic  case  with  the  same  Holder  exponent  and  degree  of  compressibility. 

In  the  present  context,  the  separation  vector  between  two  fluid  particles  is 
a  stochastic  process  solving  an  Ito  stochastic  differential  equation  with  non- 
Lipschitz  diffusion  coefficient.  To  guarantee  the  existence  and  the  uniqueness 
in  law  of  the  solution,  we  shall  add  pure  diffusion  to  the  velocity  field.  By 
considering  an  appropriate  Lyapunov  function,  we  shall  demonstrate  that 
there  exists  an  invariant  measure  for  the  fluid-particle  separation.  The  in¬ 
variant  measure  is  unique,  ergodic,  and  non-degenerate  as  a  consequence  of 
the  irreducibility  and  the  strong  Feller  property  of  the  process. 

The  paper  is  divided  as  follows.  Section [2] describes  a  generalization  of  the 
Kraichnan  model  on  a  d-dimensional  cylindrical  surface.  The  two-dimensional 
case  is  studied  in  detail  in  section  [3]  Sections  |4]  and  [5]  contain  the  results  on 
the  fluid-particle  separation  and  its  invariant  measure.  The  limit  of  vanishing 
diffusivity,  the  effect  of  a  viscous  regularization  of  the  velocity  field,  and  the 
role  of  the  Prandtl  number  are  discussed  in  section  EO 

For  d  =  4,  the  collapsed  phase  can  exist  only  for  smooth  velocity  fields  (£  =  2). 
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2  Kraichnan  model  on  a  d-dimensional  cylindrical  surface 


We  consider  the  dynamics  of  fluid  particles  in  a  turbulent  flow  on  a  d-dimen¬ 
sional  cylindrical  surface  S.  The  velocity  field  is  a  family  of  white  noises 
taking  their  values  in  the  space  of  vector  fields  on  S.  Specifically,  v(t,  x)  is  a 
Gaussian  stochastic  process  with  zero  mean  and  covariance 


E (va(t,x)vp{s,y))  =  Dap(x  -  y)S{t  -  s),  (1) 

where  x,y  £  x  [— nL,  nL)d~d'  C  (d  >  d!)  and  L  is  the  radius  of 
the  cylinder.  The  velocity  is  by  definition  statistically  homogeneous  in  space, 
stationary  in  time,  and  invariant  under  time  reversal.  Moreover,  we  assume 
periodicity  in  the  d  —  d!  “radial”  coordinates. 

It  is  convenient  to  write  the  spatial  covariance  Dap(r)  in  terms  of  its 
Fourier-space  representation: 


Dap{r) 


1  ik" -r" 

(2ir)d' (2n  L)d~d' 


dk' elk' r'  Faf}(k) 


with  k  =  (k',k")  £  R.d>  x  j-Zd~d' ,  r  =  (■ r',r ")  G  x  [—'KL1TrL)d~d', 
and  a,  /3  =  1 , ,d.  The  presence  of  a  series  in  the  fc"-coordinates  accounts 
for  the  periodicity  of  the  velocity  field  in  the  r"-coordinates. 

We  adopt  the  following  form  for  the  spectral  tensor: 


Fap{k)  = 


Aatf)  (fcj  p) 


with  £  £  R+,  £  €  [0,2],  p  £  [0, 1],  and 


Aap{k-,  p)  =  (1  -  p)dQ/3  +  (pd  -  1) 


(2) 


As  we  shall  see  in  the  latter  part  of  this  section,  Fap(k )  has  been  chosen 
in  such  a  way  that,  in  the  limits  L  — >  0  and  L  — >  oo,  Dap(r)  tends  to  the 
covariance  of  an  isotropic  random  field  @ 

The  spectral  tensor  is  real,  symmetric  (Fap(k)  =  Fpa(k))  and  non¬ 
negative  definite  Vfc  £  x  j-Zd~d  ,  i.e., 


^2  Fap(k)uaup  >  0  V(ui,...,ud)  £  Rd, 

1  <a,/3<d 


as  can  be  checked  using  the  Cauchy-Schwartz  inequality.  These  properties 
guarantee  that  Dap(r)  is  the  spatial  covariance  of  a  homogeneous  random 
field  (e.g.,  ref.  [25] .  p.  20).  Moreover,  Fap(k)  is  an  even  function  of  fc,  and 
therefore  the  velocity  is  statistically  invariant  under  parity:  Dap(—r )  = 
Dap(r). 

3  The  spectral  tensor  could  in  principle  be  multiplied  by  a  positive  coefficient 
determining  the  intensity  of  the  velocity  fluctuations.  For  the  sake  of  simplicity,  we 
set  that  coefficient  to  one. 
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As  a  consequence  of  statistical  homogeneity  and  parity  invariance,  the 
covariance  of  velocity  differences  can  be  expressed  in  terms  of  Dap(r): 

E([va(t,x  +  r)  -  va{t,  x)][vp(s,x  +  r)  -  vp(s,x)])  =  2 dap(r)5(t  -  s )  (3) 

with  dap(r)  =  Dap( 0)  -  Dap(r)  [23] . 

The  meaning  of  the  parameters  p,  £,  and  £  may  be  understood  by  con¬ 
sidering  the  limit  of  Dap[r)  for  L  — >  oo  and  for  L  — >  0. 

The  limit  L  — >  oo  (and  1/L  — >  dk")  yields: 


lim  Da0(r ) 

L—t  oo 


l 

(27r)d 


eik'-r'+ik".r"Aafj 

(Ilfc'H2 +  ||fc"||2  +  £-2)^ 


1  f  dh  eikrAap(k]fp) 

(2ttY  JRd  (||fe||2  +  ^-2)^ 


(4) 


In  this  limit,  Dap(r)  tends  to  the  spatial  covariance  of  a  d-dimensional 
isotropic  field  with  correlation  length  £  and  degree  of  compressibility  p  [15.; 
mi-  The  parameter  £/2  represents  the  inertial-range  Holder  exponent  of  the 
velocity:  X)«=i  4»(r)  =  O(IMI^)  as  \\r / 1\\  —>  0.  For  £  =  0  the  velocity  field 
is  purely  diffusive;  for  £  =  2  it  is  spatially  smooth,  and  its  spatial  regularity 
decreases  with  decreasing  £.  In  particular,  the  Kolmogorov  scaling  is  obtained 
for  £  =  4/3,  for  the  time  integral  of  <[3|)  must  be  proportional  to  |M|4 5/3  in 
Kolmogorov’s  phenomenology  (13  0 

It  is  worth  noting  that  for  a  finite  L  equation  ff}  describes  the  velocity 
covariance  at  space  separations  much  smaller  than  L. 

In  the  second  limit,  L  — >  0,  we  obtair0 


lim  Dap(r)  =  6(r") 


JC_ 

(27r)d' 


eik'  r'  Aap(k' \  p') 
(life'll2  +  £-2)^ 


with  £'  =  £  +  (d  —  d')  and 


,  =  p{d  -  1) 

^  p{d  —  d')  +  d'  —  1  ’ 


AC  =  1  + 


p(d  —  d') 
d'-l 


if  d!  >  1, 


An(fc,;p/)  =  1,  K.  =  p(d —  1)  if  d'  =  1  and  p  >  0. 

4  The  same  conclusion  can  be  reached  rigorously  by  defining  the  Kraichnan  en¬ 
semble  as  the  limit  of  an  Ornstein-Uhlenbeck  process  for  vanishing  correlation 
time  |10|. 

5  This  can  be  shown  by  multiplying  Dap(r)  by  a  function  integrating 

over  r"  £  [— ivL,nL),  taking  the  limit  L  — >  0,  and  noting  that  only  the  term 
corresponding  to  k"  =  0  has  a  non-zero  limit  equal  to 


m 


k.  r 

(2n)d'  JRd' 


dk' 


eik'r’  Aap{k'-p') 

(||fe'P  +£~2) 


d’+e’  ■ 
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We  thus  recover  the  covariance  of  a  d'-dimensional  isotropic  velocity  field 
with  Holder  exponent  £'/2,  correlation  length  S,  and  degree  of  compressibil¬ 
ity  p'  mm- 

The  exponents  £  and  £'  must  satisfy  the  inequalities  0  <  £  <  2  and 
0  <  £'  <  2.  Therefore,  the  limit  L  — >  0  makes  sense  only  in  two  cases: 

a)  d'  =  d  —  1,  £  £  [0, 1],  and  £'  =  £  +  1  £  [1,  2]; 

b)  d'  =  d  —  2,  £  =  0,  and  £'  =  2. 

We  are  interested  in  the  situation  where  the  fluid  particles  disperse  when  L  — > 
oo  (d-dimensional  isotropic  flow)  and  collapse  when  L  — >  0  ((/-dimensional 
isotropic  flow).  Moreover,  we  focus  on  spatially  rough  velocity  fields  leaving 
aside  the  cases  £  =  0  and  £'  =  2.  This  situation  can  be  realized  only  in 
case  a),  for  £  £  (0, 1),  and  under  the  conditions  [15]: 


d  ,  ,  d' 

P  <  £2  and  P  >  ^2- 

The  first  inequality  is  actually  satisfied  for  all  d  and  £  given  that  £  £  (0, 1) 
and  p  £  [0,1].  The  second  inequality  can  be  rewritten  in  terms  of  p  as  follows: 


P  > 


d-2 

w+zy 


(5) 


The  restriction  0  <  p  <  1  and  inequality  ©  imply  the  additional  condi¬ 
tion  d  <  5. 

In  the  remainder  of  the  paper,  we  shall  investigate  the  statistics  of  fluid- 
particle  separations  on  a  two-dimensional  cylindrical  surface  (d=  2). 


3  Two-dimensional  cylindrical  surface 


For  d  =  2,  case  a)  is  the  only  realizable  one,  corresponding  to  d!  =  1.  Con¬ 
dition  d5J)  reduces  to  p  >  0  independently  of  £. 

The  spatial  covariance  of  the  velocity  field  takes  the  form 


Dap{r) 


4ir2L 


OO 

E 


£);p) 


(6) 


with  r  =  (ri,  r%)  £  fl  =  Kx  [— nL,  i xL).  In  eq.  |6|  we  have  written  k  =  (fci,  ^2) 
with  ^2  =  j/L  and  j  £  Z  to  make  the  dependence  on  L  explicit.  We  shall 
keep  this  notation  in  the  remainder  of  the  paper. 

We  now  restrict  attention  to  space  separations  much  smaller  than  t.  For¬ 
mally,  this  is  equivalent  to  considering  the  limit  i  — >  00.  The  spatial  variance 
of  the  velocity  field,  Dap  (0) ,  diverges  as  l  tends  to  infinity  (appendix  [A} ; 
this  behavior  reflects  the  divergence  of  the  average  kinetic  energy  of  the 
fluid.  Nevertheless,  dap(r)  has  a  finite  limit  for  all  r,  and  the  statistics  of 
velocity  differences  remains  well  defined. 
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The  limit  of  dap(r)  for  l  — >  oo  can  be  computed  explicitly  (appendix [A}. 
The  correlation  of  the  axial  component  is  written: 


lim  dn(r)  = 


+ 


P 


r(-i±i) 


23+«7r3/2r 


(i  +  l) 


N 


i+? 


2t r3/2A  (2  +  |) 


E-? 

3= l 


/ 1 +  ^  ~  P)£  r  £ 


—  2  cos 


jr  2 


/  J>i 


P^l  + 

+(1  —  2p) 


2}  \  2L 

j  In 


2A 


r 

2 

i±i 

A'i+c 

2 

3  +  € 

A'a+f 


jN 


A 

jN 


.  (7) 


where  Kv(z)  denotes  the  modified  Bessel  function  of  the  second  kind  of  or¬ 
der  ^  and  argument  2.  The  correlation  of  the  radial  component  has  the  form: 


(1-P) 

,hm  d22(r)  = - 7 - 77— 

e^oo  23+«7r3/2r  M  +  L 

+  27r3/2r  (2  +  I)  j=i 


ni+« 

i-$  / 1  +  p£  p  1  +  £ 


—  2  cos 


JT2 

L 


(1  -  P)  1  + 


£  W ikil 


2A 


AT 


+(2p-l)(^)  A"  3+s 


3+5 


ihJ 

L 

j\n\ 


■  (8) 


Finally,  the  mixed  correlations  can  be  written  as  follows: 


lim  d\2(r)  =  lim  d2i(r) 

■C—too  £—}oo 


(2  P  —  1)A^  i_| 

27r3/2A  (2+|)  j=1 


(9) 


The  limit  £  —>  00  will  be  hereafter  understood. 


3.1  Large-scale  form  of  the  covariance  of  velocity  differences 


To  understand  the  nature  of  the  random  velocity  field,  it  is  useful  to  consider 
the  covariance  of  velocity  differences  at  space  separations  much  greater  than 
the  radius  of  the  cylinder. 

The  series  in  eqs.  ©-©  converge  uniformly  (appendix  0.  For  |ri|/L  — > 
oo,  it  is  therefore  possible  to  deduce  the  asymptotic  expansion  of  dap(r) 
from  the  limiting  behavior  of  the  single  terms  of  the  series.  The  asymptotic 
expansion  of  Ku(z)  for  z  — >  oo  is  (e.g.,  ref.  9],  formula  II  7.13(7)) 


(I  arg z\  <  3tt/2). 


Thus,  the  ^-dependent  contributions  to  dap(r )  decay  exponentially  fast  with 
increasing  space  separation.  The  remaining  contributions  give 


with 


and 


= 


In  the  latter  equation 


dn(r)  ~  £>i  |ri|1+s  +  xq  as 


2>i  = 


[i  +  (i  -  P)t)Lir  (i±£) 
271^/2(2  + or  (l  +  I) 


oo 


p 

23+?7rV2r  (l  +  1 

-)L 

C(i  +  0- 


(10) 


cw  =  E^  (s>1) 


3  = 1 


is  the  Riemann  Zeta  function.  Likewise,  we  have 


d22(r)  ~  D-2  |?"i|1+?  +  >^2  as 

±J 


with 


®2  =  — Si, 


X2  = 


1  +  p£ 


OO 


'  ^1- 


p  1  +  (i  -  p)£ 

Finally,  the  off-diagonal  terms  vanish  at  large  space  separations: 

lim  di2(r)  =  lim  c?2i(r)  =  0. 

|ri|/ L—^oo  \r\\/ L—±oo 


(ID 


(12) 


The  above  asymptotic  expressions  show  that,  at  separations  much  greater 
than  the  radius  of  the  cylinder,  the  velocity  difference  may  be  regarded  as 
the  superposition  of  two  independent  one-dimensional  random  fields.  One 
field  is  directed  along  the  axial  direction;  the  other  one  is  directed  along  the 
radial  direction.  Both  the  fields  depend  only  on  r\.  In  particular,  the  axial 
field  is  a  one-dimensional  Kraichnan  velocity  field  with  Holder  exponent  1  +  £ 
like  the  one  considered  in  ref.  [30] . 

At  large  separations,  the  small-scale  dynamics  manifests  itself  through  an 
effective  diffusivity  represented  by  the  constants  it\  and 
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4  Fluid-particle  dynamics 


In  the  present  context,  the  separation  between  two  fluid  particles  can  be 
regarded  as  a  stochastic  process  on  17  with  diffusion  coefficient  dap(r)  (and 
drift  coefficient  equal  to  zero).  To  ensure  the  (weak)  existence  and  uniqueness 
of  the  trajectories  of  the  process,  we  add  diffusion  to  the  velocity  field  and 
replace  dap{r)  by 

d%p(r)  ■=  dap(r)  +  2n5ap,  k  >  0. 

The  additional  term  can  model  the  action  of  molecular  diffusion  on  fluid 
particles  as,  e.g.,  in  ref.  m-  The  constant  k  will  be  referred  to  as  diffusivity. 

The  separation  vector  between  two  fluid  particles  will  be  denoted  by  R. 
According  to  the  above  remark,  R  satisfies  the  ltd  stochastic  differential 
equatiorO 

dR(t)  =  V2a(R(t))dB(t),  U(0)  =  r  €  17,  (13) 

where  B  is  Brownian  motion  on  17  and  a  is  defined  through  the  Cliolesky 
decomposition  of  the  matrix  dK\ 

aaT  =  dK 


with 


cr(r) 


( Vdn(r )  0  \ 


<*12  fr) 
W<*u(r) 


[dK12(r)}2  > 

<*n  (r)  / 


r  G  17. 


Although  the  diffusion  matrix  dK  is  nor  Lipschitz  continuous  nor  bounded, 
the  existence  and  uniqueness  of  the  solution  of  eq.  m  can  be  proved  using 
Stroock’s  and  Varadhan’s  theory  of  martingale  problems  [28].  To  directly 
exploit  this  theory,  we  shall  first  consider  the  periodic  extension  of  eq.  m 
on  M2,  and  then  project  the  resulting  process  on  17.  We  therefore  introduce 
the  projection  p  :  R2  — >  f2  with 


p(r)  =  (r1,-nL  +  r2-2nL  ^-J) 


and  define  a  :=  <r  op.  Likewise  we  denote  dK  :=  dK  o  p  =  ddT. 

Before  proceeding  further,  it  is  convenient  to  define  some  notation.  The 
spaces  of  bounded  measurable  and  bounded  continuous  functions  on  17  will 
be  denoted  by  J?b(l 7)  and  ^(17),  respectively.  The  set  ^2(17)  will  be  the 
space  of  functions  having  two  continuous  derivatives.  Analogous  definitions 
will  apply  to  functions  defined  on  K2. 

6  If  A  and  Y  denote  the  positions  of  two  fluid  particles,  the  separation  vector 
between  the  two  particles  is  defined  as  R  :=  Y  —  X .  The  common  physical  notation 
for  the  evolution  equation  for  R  would  be 

rJ  d 

-jj-  =  SRv(t)  +  2-v/«£(t) 

where  £  is  white  noise  and  the  statistics  of  Snv(t)  :=  v(t,Y(t))  —  v(t,X(t))  is 
defined  by  eq.  m 
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Proposition  1  The  ltd  stochastic  differential  equation  onR2: 

dR(t)  =  V2 &(R(t))dB{t),  R{ 0)  =  r  £  R2,  (14) 

where  B  is  the  standard  Brownian  motion  on  R  ,  has  a  unique  (in  law)  weak 
solution.  In  particular,  the  solution  is  a  continuous  Markov  process. 

For  r  £  R2  and  U  C  R2  measurable,  let 

P(0,  r;  s,  U)  :=  P  (r(s)  £  U  if  R(0)  =  r) 

be  the  transition  probability  distribution  of  R,  and  let  ( Tt)t>0  be  the  associate 
transition  semigroup: 

Ttf(r):=[  f(p)P(0,r;t,dp) 

Jr2 

with  f  £  Then,  the  semigroup  ( Tt)t>0  has  the  strong  Feller  property, 

i.e.  Tt(3§b{ R2))  C  %(R2)  for  all  t  >  0. 

Proof  The  diffusion  matrix  has  the  following  properties: 

1.  dK  is  continuous; 

2.  dK(r)  is  symmetric  and  strictly  positive  definite  for  all  r  £  R2.  The  dif- 
fusivity  k  is  indeed  assumed  to  be  strictly  positive,  and  the  spatial  co- 
variance  of  velocity  differences  must  be  symmetric  and  uniformly  non¬ 
negative  definite  for  all  r  (ref.  p.  97),  i.e., 

dap{r)uaup  >0  V  (u1,u2)  and  r  £  R2; 

l<a,/3<2 

3.  there  exists  a  positive  constant  C\  such  that  for  all  a,  ff  and  r 

\dap(r)\  <  Ci(l  +  ||r||2).  (15) 

This  property  is  a  consequence  of  the  asymptotic  behaviors  m  to  m 
and  of  the  the  fact  that  dK(r )  is  bounded  at  the  origin  and  continuous 
on  R2. 

Under  the  above  conditions,  Stroock’s  and  Varadhan’s  uniqueness  theorem 
apply  to  the  martingale  problem  for  dK  Then,  the  proposition  follows 
from  the  equivalence  between  the  wcll-posedness  of  martingale  problems  and 
the  existence  and  uniqueness  in  law  of  weak  solutions  of  stochastic  differential 
equations  (ref.  [27],  pp.  159  and  170).  □ 

The  process  R  can  be  regarded  as  the  projection  of  R  on  17:  R(t)  =  p  (R(t)). 
The  properties  of  R  can  then  be  deduced  from  those  of  R. 

Corollary  1  Equation  (11311  has  a  unique  (in  law)  weak  solution.  The  tran¬ 
sition  semigroup  (Tt)t> o  associated  with  R  has  the  strong  Feller  property: 
Tt(&b{H))  C  %(1 7)  for  all  t  >  0. 
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Proof  Given  an  initial  condition  r  £  17,  a  weak  solution  of  eq.  m  can  be 
constructed  by  taking  a  solution  of  eq.  m  with  the  same  initial  condition 
and  projecting  it  on  17. 

The  key  observation  to  prove  uniqueness  in  law  is  that  any  solution  of 
eq.  m  on  17  can  be  uniquely  mapped  into  a  continuous  solution  of  eq.  m 
on  R2 .  Then,  uniqueness  in  law  in  R2  guarantees  that  also  the  solution  on  17 
is  unique  in  law. 

Finally,  (Tt)t>Q  has  the  strong  Feller  property  and  the  projection  p  is 
locally  invertible  (with  continuous  inverse).  Hence,  ( Tt)t>o  has  the  strong 
Feller  property.  □ 


5  Invariant  measure  of  fluid-particle  separations 


An  invariant  measure  for  R  can  be  constructed  by  adapting  to  the  case 
under  examination  the  procedure  described  in  ref.  [261 .  Clearly,  an  invariant 
measure  may  exist  only  if  the  trajectories  of  the  stochastic  process  do  not 
“escape  to  infinity”.  To  control  the  behavior  of  the  first  component  of  R, 
which  is  not  bounded,  we  therefore  introduce  the  Lyapunov  function  V  : 
17  — »  R+: 


V(r ) 


h(h  +  1  )c4  +  2(1  —  h2)c2r\  —  h{  1  —  h)r\ 


1 


2(1  -  h) 


\n 


4(1  -  h)c2(h+1'> 
12(1— fc) 


if  |n  |  <  c, 


if  |n  |  >  c, 


where  c  >  0  and  1  >  h  >  0.  The  function  V  is  twice  continuously  differen¬ 
tiable  and  has  the  asymptotic  behavior  needed  for  the  proof. 

Lemma  1  If  (1  +£) /2  >  h  >  1/2  and  A  denotes  the  infinitesimal  generator 
of  (Tt)t> 0,  then  the  Lyapunov  function  has  the  following  properties  for  all  £  £ 
(0,1):" 

1.  lini||I,||_>00  AV(r)  =  -oo; 

2.  there  exists  m  €  R  such  that  AV{r)  <  m  for  all  r  £  17; 

3.  TtV(r)  =  V(r )  +  [  TsAV(r)ds. 

Jo 

Proof  For  /  £  ^2(17),  the  infinitesimal  generator  of  (Tt)t>0  has  the  form 


Af(r)  =  tr  [cr(r  )erT(r  )D2  f(r )], 

where  D2  f  denotes  the  Hessian  of  the  function  /.  The  action  of  the  genera¬ 
tor  A  on  V(r)  is  written: 

AV(r)  =  a211(r)^-. 

From  eq.  m  we  obtain 


AV(r)  ~  -(2 h  -  l)Si|n|1+5_2'1 


as  || r- 1|  — >  oo. 
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If  (l  +  £)/2  >  h  >  1/2,  we  have  limn^n^oo  AV(r)  =  —  oo  for  all  £  €  (0, 1).  It  is 
worth  noting  that  this  latter  result  relies  on  the  fact  that  dn(r)  =  0(|ri|1+^) 
as  ||r||  — >  oo  with  1  +  £  >1. 

Property  2  is  a  consequence  of  the  continuity  of  AV  and  of  property  1. 
Finally,  the  transition  semigroup  satisfies 

Ttf(r)  =  f(r)+  f  TsAf(r)ds 
Jo 

for  any  /  £  ^b2(  12 ).  The  same  property  holds  true  for  the  function  V,  as  can 
be  shown  using  an  approximation  procedure  similar  to  the  one  described  in 
ref.  [2BJ,  pp.  167  -168.  The  details  are  given  in  appendix  iBl  □ 

We  now  make  use  of  the  properties  of  the  Lyapunov  function  to  obtain  the 
following  result. 

Proposition  2  There  exists  an  invariant  measure  p  for  the  stochastic  pro¬ 
cess  R,  i.e. 

[  {Ttf)(r)n(dr)=  [  f(r)p{dr)  (16) 

J  Q  J  Q 

for  all  f  £  and  for  all  t  >  0. 

Proof  Given  ro  £  12,  we  have 


1 

t 


AV)(r0)ds 


TsAV(r0)ds 
|  V(r0)-TtV(r0) 


<m  + 


V(ro) 

t 


Hence 


We  now  introduce 
for  any  /  G 


sup 
t>  l 


1 

-  /  Ts(m-  AV)(r0)ds 
t  J  o 


<  oo. 


the  family  of  “average”  measures  (pt)t> l 
as 


(17) 


on  12  defined, 


[  f(r)nt(dr)  =  7  [  Tsf(r0)ds. 

Jo  1  Jo 

We  show  that  (/q)t> i  is  uniformly  tight.  For  a  given  N  >  0  we  define 
En  =  {r  £  12  :  to  —  AV(r)  <  N}. 


The  set  En  is  compact:  it  is  closed  since  it  is  the  preimage  of  a  closed 
subset  of  12,  and  must  be  bounded  since  AV(r)  — >  —  oo  as  ||r||  — >  oo.  As  a 
consequence  of  Markov’s  inequality,  the  measure  of  the  complement  of  En 
satisfies: 


AV(r)]fj,t(dr) 


1 

m 


Ts(m 


AV)(r0)ds. 
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Using  eq.  m  we  conclude  that  for  all  e  >  0  there  exists  En  C  fl  with 

i  r  i  /■* 

N=  —  sup  -  /  Ts{m  —  AV)(r0)ds 
e  t>i  [t  Jo 

such  that  ptfEf,)  <  e  for  all  t  >  1.  The  family  (pt)t>i  is  therefore  uni¬ 
formly  tight.  Then,  there  exists  a  measure  p  and  a  sequence  (tp)p>o  with 
linip^oo  tp  =  oo  such  that  ptp  converges  weakly  to  /i  as  p  — >  oo.  This  means 
that  jQfdptp  — >  f  dp  as  p  — >  oo  for  all  /  G  %,{Q)  (e.g.,  ref.  [6],  theo¬ 
rem  11.5.4,  p.  404). 

We  now  show  that  p  is  invariant.  For  any  /  G  ^,(1?),  for  t  >  0,  and  for 
all  p  such  that  tp  >  t,  we  have: 


m 


f(,r)ntp(dr)-  [  Ttf(r)ptp(dr) 
J  n 


1 


-  Tsf(r0)ds~-  TsTtf{r0)ds 


tp  Jo 

ir 

t. 


Tsf(r0)  ds  - 


Lp  JO 

1 


f't+tp 


P  Jo 


Tsf(r0)ds 


Tsf(r0)ds  - 


P  Jt 

1 


» t-\-tp 


P  JO 


Tsf(r0)  ds 


P  Jtp 


< 


2*||/||c 


Hence,  for  all  t  >  0, 


[  f(r)Htp(dr)~  [  Ttf(r)  fitp(dr)  ^  0  as  p  ->•  oo. 

J  n  J  n 

The  semigroup  ( Tt)t>o  satisfies  Tt(f&b(fi))  C  %,{Q)  since  (T)t> 0  has  the 
strong  Feller  property.  By  using  the  weak  convergence  of  ptp  to  p,  we  can 
thus  conclude  that  (fTCl)  holds  for  all  /  G  %,{Q)  and  for  all  t  >  0.  The 
measure  p  is  therefore  invariant  for  R.  □ 

To  show  that  the  invariant  measure  is  actually  unique,  we  need  the  following 
result  stating  that  R  has  no  closed  invariant  set  different  from  the  whole 
space. 

Lemma  2  The  semigroup  ( Tt)t>o  is  irreducible,  i.e.  the  transition  probabil¬ 
ities  of  R,  P( 0,  r;  t ,  U),  are  strictly  positive  for  all  t  >  0,  for  all  r  G  fi,  and 
for  all  non-empty  open  sets  U  C  17. 

Proof  For  all  r  G  R2  the  linear  transformation  associated  with  a(r)  is  in¬ 
vertible,  and  therefore  maps  R2  into  itself.  Hence,  the  semigroup  (Tt)t>0  is 
irreducible  (ref.  [29],  theorem  24,  p.  66). 

The  transition  probabilities  of  R  are  connected  to  those  of  R  as  fol¬ 
lows:  P(0,r;t,U)  =  P(0,r*\t,U*)  where  U*  =  p  1  (C/)  and  r*  is  any  point 
in  p~1({r}).  Therefore,  {Tt)t>o  is  irreducible.  □ 

We  can  now  state  the  main  result  regarding  the  invariant  measure  of  R. 
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Theorem  1  There  exists  a  unique  invariant  measure  p  for  the  stochastic 
process  R.  The  measure  p  is  ergodic  and  equivalent  to  any  transition  proba¬ 
bility  P(0,r;t,U)  with  r  £  R,  t  >  0,  and  U  C  17  measurable.  Moreover,  p  is 
absolutely  continuous  with  respect  to  the  Lebesgue  measure,  and  is  therefore 
non  degenerate  (i.e.,  broad  in  r). 

Proof  We  have  already  proved  that  p  is  invariant.  Its  uniqueness,  ergodicity, 
and  equivalence  to  any  transition  probability  follow  from  the  fact  that  the 
transition  semigroup  associated  with  R  has  the  strong  Feller  property  and  is 
irreducible  Bin]  (see  also  ref.  j5] ,  chapter  4) . 

To  prove  the  absolute  continuity  of  p  with  respect  to  the  Lebesgue  mea¬ 
sure,  we  introduce  the  family  of  transition  probabilities 

pOO 

Qx{r ,  U)  =  X  e~XsP( 0,  r;  s,  U)  ds 

Jo 

with  r  G  i?  and  U  C  i?  measurable,  as  well  as  the  associate  transition 
semigroup 

Tx f(r)=  f  f(y)Q\(r ,  dy). 

Jn 

Likewise,  we  define  an  analogous  family  Q\(r,U)  for  the  process  R.  The 
measure  p  is  invariant  also  for  (7a)a>o: 

[  p(dy)T\f{y)  =  A  [  dse~Xs  [  p{dy)Tsf(y)  =  [  p{dy)f(y) 

J  £2  Jo  J  Cl  J Cl 

for  any  /  €  SSh{TT),  and  hence 

h(U)=[  p(dy)Q\(y,  U)  (18) 

J  a 

for  any  measurable  set  t/  C  17. 

For  all  r  £  R2,  the  measure  Q\(r,-)  is  absolutely  continuous  with  re¬ 
spect  to  the  Lebesgue  measure  (see  ref.  [25],  theorem  10,  p.  24).  It  fol¬ 
lows  that  Q\{r,-)  has  the  same  property  for  all  r  £  17  since  Q\(r,U)  = 
Qx(r* for  a  given  r*  £  p-1{r}.  From  eq.  (1151).  p  is  therefore  abso¬ 
lutely  continuous  with  respect  to  the  Lebesgue  measure.  □ 

As  a  consequence  of  lemma  [2]  and  theorem  [TJ  the  transition  probability  of  R 
has  a  positive  density  with  respect  to  the  Lebesgue  measure:  P( 0,  r;  t,  dp)  = 
p(0,  r;  t,  p)dp.  The  probability  density  function  is  the  (possibly  weak)  solu¬ 
tion  of: 

dtp  =  Mp ,  (19) 

where,  for  f  £ff2(i 7), 

MCP)=  dp*dPpd«p{p)f(p)- 

l<«,/3<2 
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6  Conclusions 


We  have  studied  the  dynamics  of  fluid  particles  in  a  compressible  turbulent 
velocity  field  on  a  cylinder.  The  model  that  we  have  introduced  is  a  general¬ 
ization  of  the  isotropic  Kraichnan  ensemble.  Although  the  parameters  of  the 
velocity  have  been  set  in  such  a  way  as  to  produce  explosive  separation  of  the 
fluid  particles  in  the  isotropic  limit  (L  — >  oo),  on  the  cylinder  the  probability 
distribution  of  the  separation  tends  to  an  invariant  measure.  This  behav¬ 
ior  is  a  result  of  the  compressibility  effects  generated  at  large  scales  by  the 
compactification  of  the  “radial”  dimension. 

The  diffusivity  n  has  been  taken  strictly  positive  to  guarantee  the  exis¬ 
tence  of  solutions  to  eq.  m-  The  addition  of  Brownian  motion  to  Lagrangian 
trajectories  influences  the  dynamics  of  fluid  particles  at  small  separations. 
Therefore,  the  presence  of  a  nonzero  diffusivity  may  be  relevant  for  the  non¬ 
degeneracy  of  the  invariant  measure,  but  should  not  affect  the  existence  of 
the  invariant  measure  itself,  which  rather  depends  on  the  large-scale  form 
of  the  velocity  field.  The  limit  k  — >  0  may  be  tackled  by  means  of  Wiener 
chaos  decomposition  methods  We  conjecture  that  our  results 

remain  valid  in  that  limit.  Indeed,  in  the  situation  considered,  the  small-scale 
dynamics  of  fluid  particles  is  the  same  as  in  the  weakly  compressible  phase 
of  the  isotropic  Kraichnan  ensemble.  In  that  regime,  Lagrangian  trajectories 
separate  in  time  even  for  vanishing  k  owing  to  the  poor  spatial  regularity  of 
the  velocity  [15].  Thus,  the  invariant  measure  should  remain  non-degenerate 
as  k  — ^  0. 

The  Reynolds  number  is  infinite  in  our  study  since  the  viscosity  of  the 
fluid,  v,  has  been  set  to  zero  from  the  beginning.  For  the  same  reason  the 
Prandtl  number  Pr  =  v/n  is  equal  to  zero.  The  viscosity  can  be  taken 
into  account  by  multiplying  the  spectral  tensor  (O  by  the  factor  e~v  , 
where  77  oc  v3!4  plays  the  role  of  the  viscous  length  of  the  flow  m-  This 
modification  has  a  small-scale  regularizing  effect  on  the  velocity  field,  which 
for  any  77  >  0  is  locally  Lipschitz  continuous.  Obviously,  a  positive  77  does  not 
alter  the  proofs  of  the  results  shown  in  the  paper. 

The  order  of  the  limits  n  — >  0  and  77  — >  0,  however,  deserves  a  detailed 
discussion.  Taking  the  limit  k  — >  0  before  77  — >-  0  is  equivalent  to  letting  Pr 
tend  to  infinity.  The  opposite  order  corresponds  to  the  limit  Pr  — >  0.  As 
first  observed  in  ref.  |8],  when  these  limits  are  considered  the  range  of  weak 
compressibility  splits  into  two  ranges:  what  is  now  called  the  range  of  weak 
compressibility  in  the  strict  sense,  0  <  p  <  {d—  2  +  £)/  (2£),  and  the  range  of 
intermediate  compressibility,  (d— 2+£)/(2£)  <  p  <  d/£2.  In  the  former  range, 
the  order  of  the  limits  n  — >  0  and  77  — >■  0  is  not  relevant  for  the  Lagrangian 
dynamics  [II].  At  small  scales,  fluid  particles  disperse  irrespective  of  the  order 
of  the  limits,  and  therefore  we  expect  the  invariant  measure  of  the  separation 
to  be  non-degenerate.  By  contrast,  the  order  matters  in  the  latter  range  [14|. 
For  intermediate  values  of  the  compressibility,  if  the  viscous  regularization  is 
removed  before  the  diffusivity  (Pr  — >  0),  the  small-scale  Lagrangian  dynamics 
is  once  more  characterized  by  the  explosive  separation  of  the  trajectories.  If  n 
goes  to  zero  before  77  (Pr  — »•  00),  the  trajectories  coalesce  also  at  small  scales, 
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and  the  invariant  measure  of  the  separation  should  degenerate  into  a  Dirac 
delta  function. 

In  summary,  we  believe  that  the  present  study  captures  the  behavior  of 
the  Lagrangian  trajectories  on  cylindrical  manifolds  for  all  Pr  except  for  the 
limit  Pr  — >■  oo  in  the  intermediate-compressibility  regime.  These  results  are, 
moreover,  relevant  to  turbulent  transport  of  passive  scalar  fields  in  virtue  of 
the  relation  subsisting  between  the  scalar  correlations  and  the  dynamics  of 
fluid  particles  H5KU- 

We  conclude  by  noting  that  when  both  the  dimensions  of  the  plane  are 
compactified  one  obtains  the  Kraichnan  flow  on  a  two-dimensional  periodic 
square  studied  in  ref.  [4].  The  velocity  field  considered  there  was  however 
smooth  in  space. 


Acknowledgements  The  authors  are  grateful  to  F.  Flandoli  and  Y.  Le  Jan  for 
fruitful  discussions. 


A  Covariance  of  velocity  differences 


For  d  =  2,  d!  =  1,  £  £  (0, 1),  £  (1,  2),  the  spatial  covariance  of  the  velocity  field 

takes  the  form 


Dap(r ) 


1 

47r2L 


e<fcirMa/3((fc1,f);p) 

(fci  +  jj  +  ji) 


with  reft  We  first  establish  the  convergence  of  the  above  series.  It  is  convenient 
to  denote  the  integrals  by  D^(ri)  and  thus  rewrite  the  covariance  as  follows: 


Dap{r ) 


1 

4t PL 


E  D^yir2 


j  =  —  oo 


Using  the  inequality 


<Ei) 

47r2L 


2tt2L 


cos 


3=1 


(20) 


\Aap{k)\  p)\  <  l-p  +  \2p—  1| 


Vq,/3  =  1,2  and  VfeG' 


we  obtain  that  the  coefficients  of  the  series  satisfy  for  all  r\ 
|D^(n)|<(l-p+|2p-l|)M, 
with  (e.g.,  ref.  [9],  formula  I  1.5(2)) 


Mj  = 


(i , 

r(i  +  |)  \l2  e) 


The  series  X)yLi  Mi  converges  for  all  £  £  (0, 1)  and  l  >  0  as  well  as  in  the  limit  l  — > 
oo,  as  it  can  be  checked  by  means  of  the  integral  test.  Then,  the  Weierstrass  criterion 
guarantees  that  the  series  in  the  right-hand-side  of  eq.  cap  converge  uniformly  and 
absolutely  on  ft.  The  uniform  convergence  will  allow  us  to  compute  lim^oo  dap(r) 
by  exchanging  limit  and  summation. 
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The  basic  analytical  ingredient  to  derive  eqs.  0  ®  is 


/ 


(fc2+22)-+l/2 


■dk\  = 


21-1V/2AT(|zr1[)  [nr 
r{v  +  \)\z\“ 


(21) 


with  Re{v )  >  —1/2  and  |  argz|  <  7r/2  (e.g.,  ref.  [5],  formula  II  7.12(27)). 

We  first  compute  the  correlation  of  the  axial  component  of  the  velocity;  the 
correlation  of  the  other  components  may  be  easily  derived  from  Dn(r). 

In  the  limit  t  — >  oo,  we  have 


lim  £>^0)(ri)  = 


f 


dki 


eikiri  [(l-p)£  +  Pkl] 


a  -  2p)  £■  / 


dfci 


„ik1r1 


(*?+£) 5 


+  p  j  dk\ 


k?eikiri 


+  P 


f 


dfci 


ifcir-i 


(fcl  +  1^) 


i+4 


=  (1  -  2p) 


iii7r1/2|ri|(3+«)/2A^(|iri|) 


+  p 


2(i+«)/2  r  (2  +  |) 

_i+i  2(1-«/V/2|r1|(1+«/2A-1±i  (|in|) 


^(1  +  1) 

To  compute  (0)  we  can  use  the  asymptotic  expansion  of  A+(z)  for  z  — »  0 

K+>~P(fr  +  P^(fP°+"'>  +<D  (22) 


An)~^(fr-^(ff"+o(n  d<-<2).  (23) 

Hence  we  obtain 

•  -l-g  l/2r  /1+Q 

lim  D[[* °)(0)  =  p  3  ^  2  ’ 


r{  1  +  1) 


L 

-!-e 


+  (1  -  2p) 


-i-«  7t1/2A(2±£) 

r(2  +  |) 

2A  (2+|) 


For  /  =  0  and  £  <  00,  we  have 


,(o) 


P)  =  P  J 


dk\ 


ik1r1 


-00  (A:f+£~2)22€  r(l  +  |) 


2  2  7T1' 2  ,  ,  i±i 


I3ni±+  n+s  (|(f|)- 


The  asymptotic  expansion  (f22l)  shows  that  £>^(0)  diverges  like  £1+^  as  £  — >  oo.  By 
using  expansion  J21,  it  is  nonetheless  possible  to  show  that 


_.l/2  I  V  f-i+il 

lim  [^^(O)  -  d[°}  (n)]  =  p - ! — ) — 

£—> oo  1  11  21+«r(i  + 1) 


2  /i  ,riii+e 
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(note  that  F(— (1  +  £)/2)  <  0  for  0  <  £  <  1).  Hence,  the  covariance  of  the  axial 
component  of  the  velocity  difference,  dn(r),  has  a  finite  limit  as  i  — >  oo. 

For  the  other  components  we  have: 


D 


U) 

22 


dk\ 


eikiri  [(1  -p)kj  +  p£\ 

2  +  g 

(kj  +  fr)  (k\  +  jj  +  -p) 


and 


D 


U) 


r° 

(ri)  =  D$(n)  =  ( 2p-l)£  J 


dki 


kie%‘ 


(k\  +  fr)  [k\  +  fr  +  jt) 


Therefore,  (n)  can  be  derived  from  _D^;(n)  by  replacing  p  with  1  —  p: 


U)  i 


lim  D^°\r i)  =  (2p  -  1) 

■t—too 


+  (!  -  P) 


(Ifni)  |n|(3+e)/2 
2(i+0/2r(2+  f) 

i±i  2(1-«)/27r1/2Xi±i(|fn|)  |n|(1+5)/2 

^(1  +  1)  : 


lim  O)(0)  = 


-i-£ 


2F  (2+|) 


lim  [D^}(0)  -  T>22'(ri)]  =  (1  ~  P) 


(0)/ 


r!/2  |_p  (■_!+£ 


21+er  (1  +  f) 


2  )  I  ,ri|i+£. 


The  mixed  correlation  can  be  obtained,  for  j  f  0,  by  differentiating  formula  m 
with  respect  to  n  and  by  using  [xv Kv(x)\  =  —xvKv- i(x)  (e.g.,  ref  [9],  for¬ 
mula  II  7.11(21)): 


lim  D 

t — >-00 


Ufto) 

12 


(n) 


i(2p-l)^ 


L 


i±i 

2 


^|n|(3+0/2  „ 

— TT7 -  A  1  +  g 

2^r(2  +  |)  ^ 


(Ifni)  sgn(n). 


Hence 


For  j  =  0  we  have 


lim  D^0)  (0)  =  0. 
.£—>•00 


L>)3)(ri)  =  0  Vn  el. 


Finally,  eqs.  0  ©  may  be  derived  by  recalling  that 

1  " 


dMr)  =  E  -  D(>i)eir2\ 


and  using  the  uniform  convergence  of  the  series  in  the  right-hand-side  of  eq.  (|20j>. 
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B  Proof  of  Lemma  Q] 

To  prove  property  3  of  lemma  [I]  we  first  observe  that 

E  (||R(f)||2)  =  ||R(0)||2  +  2E  ||cr(J?(s))||2ds^  <C2(l  +  J*  E  (||i?(s)||2)  ds 

where  C2  >  0  and  ||cr||  :=  [Tr(acrT)]  1'/“.  The  above  inequality  is  a  consequence 
of  (1151) .  Gronwall’s  inequality  then  yields 

E(||-R(t)||2)  <  C2e°2t  (24) 


for  all  t  >  0. 

Following  ref.  [25],  we  consider,  for  all  7  >  0,  the  function  :  R+  — >  R+  with 


V  7O 


-/ 


0  <  2  <  7 


)  =  J  * 

\<P7(7  +  1)  7  +  1  <z 


and  tp-,  G  l^’2(R+)  and  monotonically  non-decreasing.  Moreover,  we  define  Vy  := 
1 p-y  o  V.  Applying  A  to  Py  and  taking  into  account  (1 101)  yield 


|7lF7(r)|  =  an(r) 


^(nr))^+^(P(r))(^)2  =o(|r1|3+^)  (25) 


with  3  +  £  —  4h  <  2  as  ||r||  — >  00. 

Since  V-y  G  i#”2  (O),  we  have  for  all  7  >  0 


TtV-y{r )  =  V7(r)  +  [  Ts(AVy)(r)ds. 

Jo 

We  now  show  that  each  term  of  the  above  equation  tends  as  7  — »  00  to  the  corre¬ 
sponding  term  in  property  3  of  lemma  [T| 

Obviously,  lim7_>00  V-y (r)  =  V(r)  for  all  r  G  12.  Likewise,  V-y(R(t. ))  /*•  V(R(t)) 
almost  everywhere  as  7  — >  00.  Therefore,  by  the  monotone  convergence  theorem 
lim7_).00  TtVy(r)  =  TtV(r).  Finally,  lim.y_i.oo  AV~,(R(t))  =  AV(R(t))  and,  from 
eqs.  (1251)  and  (1221).  |AV7(.R(t))|  <  Cs(l  +  ||i2(£)||2)  with  C3  >  0  and  E  (||f?(f)||2)  < 
00.  Then,  it  follows  from  the  bounded  convergence  theorem  that 

lim  [  Ta{AV-y){r)ds  =  f  Ts(AV)(r)ds. 

7->°°  Jo  Jo 

This  concludes  the  proof. 
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